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ABSTRACT 


A preblem of routing a U.S. Coast Guard buoytender to 
service aids to navigation is formulated asa Symmetric 
traveling salesman prcblem. A heuristic algoritom is devel- 
oped which seeks the minimum distance tour which can be 
taken by the buoytender to visit the aids to navigation. A 
user's guide is provided. 

The algorithm is programmed in Convergent Iechnologies 
FORTRAN for use onthe Coast Guard Standard Terminal. 
Several proklems are solved by the algorithm producing solu- 
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I. INTRODUCTION 

Cne of the U.S. Coast Guard's missions is maintaining 
aids to navigation in the waters of the United States. 
These aids reguire periodic servicing to ensure they are on 
station and they are showing their proper watch characteris- 
tics. The responsibility of maintaining many of these aids 
falls on the Coast Guard's fleet of buoytenders. A buoy- 
tender in its area of responsibility may have as many as 200 
aids which it is required to maintain. 

The buoytender has scheduled ATON fcuns (Aids TO 
Navigation) and emergency outages it must handle duringa 
fiscal year. An emergency outage occurs when an aid is 

reported by a mariner to be showing improper characteris- 

tics or is missing from its’ station. These outages are 
handled shortly after they are reported. Several times 
during the year the buoytender plans an ATON run to handle 
the periodic servicing of aids, to ensure that they maintain 
their preper watch characteristics. 

Prior to a scheduled run the tender will r2ceive SANDS 
forms from the unit's district office for each of the aids 
which are due for service or relief. SANDS is a database 
where information is kept on a district's navigational aids. 
The district's aids to navigation branch identifies aids 
which are due for service or relief and notifies the tender. 
The district provides a list of aids which need to be 
serviced during the scheduled run and of work needed on each 
ue Ge The tender then makes plans accordingly to schedule 
and complete the work. 

The number of aids serviced on a particular scheduled 
ATON run may be as few as 5 (a day of local work) to aS many 


as 70 (an extended four-week trip servicing the west coast 


of Alaska from the Aleutian IsSlands to Point Hope). The 
Maintenance of the floating aids and fixei structures 
entails the travel of the tender from aid to aid. The tend- 
er's Cperations Officer usually plans the initial route to 
be travelled to service the aids due for routine maintenance 
and presents his profesal to the Commanding Officer, who may 
then modify this projosal. 

There are several factors which go into the decision 
process to produce tke desired route. Two of these factors 
are (1) which aids are scheduled for service, and (2) the 
distances between the aids. The desired route is usually 
the route of shortest distance which visits all the sched- 
uled aids and returns the tender to her point of origin. 
This proklem is a classic problem which operations analysts 
normally call ‘the traveling salesman problem.' Present 
route selection methods are (1) traditional or previously 
followed routes and (2) routes created by Operations 
Department personnel sitting down with chart and rule and 
selecting a route. 

It is the intent of this thesis to propose and igrplement 
a computer assisted approach as an alternative method for 
the kuoytender route selection problen. The succeeding 
Chapters of this paper contain a discussion of the trav- 
eling salesman problem and how it relates to th2 buoytender 
problem, a brief survey of possible solution methods to the 
problen, adescription of the solution process of the 
selected method, and the results of the route selection 


Program with conclusicns. 


II. NATURE OF THE PROBLEM 

In this chapter the principles of the traveling salesman 
problem will be presented and it will be shown how the buoy- 
tender problem can te described as a traveling salesman 
problem. The discussion will include a general description, 
definitions about solution methods, and the solution 
requirements of the problen. Further, a description of the 
computing hardware on which the selected solution method is 
programmed is provided. This provides an idea of the cafa- 
rilities of the hardware since the hardware has a bearing on 


the selection of the solution method. 


A. DESCRIPTION OF THE TRAVELING SALESMAN PROBLES 


A brief description of the traveling salesman problen 
(or TSP) is as follows. A salesman has n cities he must 
visit. The salesman starts in one of the cities and must 
travel so that he enters and departs each of the other 
cities only once. Upon visiting the last city the salesman 
will return to the city from which he originally started his 
travels. The desired solution to the problem is the route 
the salesman can traverse which is of the shortest total 
distance possible. 


It is usually assumed that a TSP can be associated with 


some n by n distance matrix. The elements of the matrix, 
dij. aeemeGastances LEQM cary 2 to City j, {1=1,-..-,n), 
(J=1,---,nh) ,and where djjzta, Daly scenes The traveling 


salesman is not allowed to leave city i andto return to 
Gicy 1 if his tour. 
Rh TSF can be associated with a graph, also. Ny (eplagn'g 


which the saleSman must visit is represented by a node in 
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this grarpa. There is an arc (1,3) in the graph wath length 
dj; for each node. 

A distance matrix may be classified aS sparse or dense. 
A sparse matrix will have a majority of its entries equal to 
Pht PM e ye The corresponding graph 1S sparse since it has 
few edges. A dense matrix is a matrix where a majority if 
not all its entries, with the exception of the d¢g‘'ts, are 
less than infinity. The number of possible permutations of 
Cities for a dense matrix is of order (n-1)!. Computation 
time and memory required for solution by total enumeration 
grows rapidly with the size of n due to the large number of 
possible solutions, each permutation of cities being a 
feasible tour and a possible solution. 

The distance matrix used by the buoytender for route 
selection will be a dense matrix. With the large number of 
permutations or possible solutions to the problen, sone 
means other than total enumeration must be used to solve the 
rkuoytender problem. 

The traveling salesman problem nay be broken intc two 
classés, symmetric and asymmetric. In a symmetric TSP the 
dij =4ji for all ifj, while in an asymmetric TSP ethiswaee 
not ke true. Further, a TSP may or may not be required to 
satisfy the triangle ineguality. Consider three cities i, 
J, Ke The triangle inequality states that 

dij < dik + dkj for all if#j#k#i. 
The conditions which make a TSP asymmetric, the unequal 
distances between cities, may cause the problem to violate 
the triangle inequality. 

The ISP, both symmetric and asymmetric, belongs to the 
NP Complete class of problems [Ref. 1,2]. This has one 
Major implication with respect to the complexity of solving 
the ISP; no polynomial-time algorithms are known or seen 
likely tc be devised for exact solution of the TSP. Only 


exponential-time exact algorithms are known. 
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There have been many algorithms proposed for the solu- 
tion of the TSP, each algorithm showing some advantages and 
disadvantages to its method of attacking and solving the 
problen. Several of the published algorithms are tailored 
to a farticular size of TSP. The TSP may be categorized by 
size aS either small or large. A reference could not be 
found in the literature which would give a hard numerical 
value by which the size of a TSP could be judged, but a 
small TSF is generally considered to be less then 15 cities. 
A large ISP is often considered to be in the neighborhood of 
45 or more cities. There are solution methods which solve 
the TSP exactly for small up to large problems. There are 
some methods which can solve some large problems exactly, 
but as a rule large problems are solved with heuristics or 
approximate methods. 

The circumstances of the buoytender problem satisfy the 
conditions of aclass of the traveling salesman  froblen. 
The Fuoytender problem can be characterized as a symmetric 
ISP, i.e., the distance traveled between any two aids is the 
same regardless cf direction travelled by the tender. The 
size of Euoytender problem will vary with the number of aids 
scheduled to be serviced on a specific ATON run; therefore 
the solution method for the buoytender problem will have to 
handle beth small and large problems. 

The triangle inequality as applied to the buoytender 
problem holds true, Eut how it holds needs some 2xplanation. 
With all distances positive, the triangle ineguality states 
that given three points A, B, and C, the distance from A to 
C must ke less than or equal to the distance from A to B 
plus the distance from B to C (see Figure 2.1). In the 
buoytender's problem the straight-line distance from A toC 
may bre less than (A to B) + (B to C), but some situations 
will reguire that A to C egual (A to B) + (B to C) (see 
ErgucEee 2.1), 


dac = dag + dac ° 
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Figure 2.1 Triangle Ineguality / Obstructed A to C Path 


The buoytender may be unable to take the direct coute from A 
to C due to shoal water or an obstruction. To get from A to 
C the tender will have to take the routes A to B and B to C. 
This is not an unusual Situation for buoytenders since the 


aids being serviced cften mark shoal water or obstructions. 


Be DESCRIPTION OF THE SOLUTICN REQUIREMENTS 


There are many proposed solution methods to the TSP and 
they may be divided into two groups: exact metheds9 and 
approximate or heuristic methods. An exact method or algo- 
rithm has the property that upon termination of the algo- 
rithno, the user will have the best possible or global 
optimuz solution. Heuristics or approximate methods, on the 
other hand, terminate with a local optimum solution or the 
best solution found thus far. This local optimum solution 
may be the global optimum but usually there is no proof that 
the global optimum has been found. Heuristics OF apErexci. 
mate methods solve Ey checking a subset of the numerous 
PosSible solutions and terminate when the subset has kLeen 


searched for the best answer it contains, a specified time 


ie 


limit for the run has been reached, or the memory capacity 
of the computing machine has been reached. 

A further consideration in selecting a method to solve 
the TSP is that there should be very little interactive work 
for the user. The user interactive work necessary to 
operate the program should be understandable to a majority 
of the intended users, who in this case are the buoytender's 
Operations aided Wersja or Navigator, and Quartermaster 
personnel. By requiring the interactive work to be mininal 
and simple in nature the possiblity of operator introduced 
error will te reduced. Further, if the interactive work is 
of a simple nature, the program will be useable by a 
Majority of the intended users without extensive training. 

Although it is desired to obtain the global optimum, it 
is not mandatory. The computer route selection method is to 
be a tool used in planning, andthe exact solution is not 
critical tothe operation of the buoytender. By using 
heuristics, computer time and memory may be reduced since a 
heuristic generally searches only a subset of all possible 
solutions. If a computer solution can be produced which is 
Fetter than the traditional route or a route produced by 
hand calculation, then a goal has been met, that goal being 
reducing the total distance traveled by the tender thereby 
Saving fuel and time. The computer solution 1s then a 
viable alternative tc be added to the present methods of 


route selection. 


C. CCMPUTING HARDWARE 


The computer on which the selected solution method is 
programmed is the Coast Guard's C3 Standard Terminal. The 
reascn for selecting this computer is that it is presently 
keing installed throughout the Coast Guard as a primary 
unit-level computer. Most intend users should have access 


to a C3 Terminal. 
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The C3 Standard Terminal is a standalone 16 bit micro- 
computer. The terminal comes in two basic configurations, 
the Integrated Work Station (IWS) and the Application Work 
Staticn (AWS). The IWS iS a master or standalone unit 
eguipped with a hard disk and 8 inch floppy mass storage 
system, a 390/1200 baud modem, a printer, and one of three 
RAM memory configurations - 256k, 384k, or 512k. The IWS 
may ke programmmed in BASIC, FORTRAN, COBOL, Pascal, and 


Assemrly. The AWS is more limited in its capabilities than 
the IWS. The AWS ais limited to 256k RAM and must be 
networked to a IWS station to access the peripherals. With 


a program in RAM the AWS will operate asa standalone 
computer, but does not have the power or computational speed 
of the IS. 

The solution method selected for programming will have 
to perform within the limitations of the C3's memory. The 
C3's RAM memory will be a critical factor in selection of 
the solution method. Since the computer 1S a purchased 
standalone and not a leased system, CPU time should not bea 
Critical taetor. It 1S still desired, though, to keep CPU 
time from krecoming excessive since the computer is also 
needed fcr cther work. The computer could be allowed to run 
overnight (off duty hours) for large problems, and therefore 
a soluticn method reguiring more than 12 hours to provide an 
acceptable solution would be excessive. A solution method 
which provides an answer in under an hour could be run 
during duty hours with very little impact on other users, 
particularly if it was run on a slave terminal. 

In the next chapter we take a brief look at some 
suggested sclution methods for the TSP and their applica- 
bility to the buoytender problen. 
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III. A BRIEF SURVEY OF SOLUTION METHODS FOR FHE ISP 

In this chapter a brief survey of the general solution 
methods for the TSP and some of thelr associated algorithms 
will be presented. The number of algorithms proposed for 
solving the TSP is extensive. To present all of the avail- 
able algorithms or to discuss them in detail is worthy of a 
dissertation in itself. This survey of solution methods is 
made Fecause a combination of a couple of the methods will 
ke used to solve the buoytender problen. This survey will 
ke brief in nature presenting the general solution methods 
and briefly describing some of the more popular of published 
algorithms which fall under these solution method headings.!} 

The solution methods may be classified under four 
general headings. These headings are Tour Building, Subtour 
Elimination, The One-Tree Formulation, and Tour Improvement. 
These headings describe how the traveling salesman problen 


1S apfyroached. 


A. TCUR BUILDING 


Tcur Building solution methods to construct a tour using 
the n nodes and the available arcs in the problem. Tour 
Building contains both exact and approximate algorithms. 
Some of the approaches taken by Tour Building Algorithms to 
solve the TSP are dynamic programming, ‘branch and bound' 


methods, and tour construction heuristics. 


Sees eae ae at oe 2S 2s =) 2s ees ee =e ae 


lThree excellent articles [Ref. 1,3,4] have been 
publisked which describe the TSP, ‘discuss general solution 
methods and review many of the algorithms which sclve it. 
These three articles review the reSearch conducted on the 
Seer souethe mid 1950's up to 1983. 
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Held and Karp [Ref. 5] present a dynamic programming 
approach which is suitable for small problems and can be 
modified to an approximate method for larger proplems. The 
dynamic programming approach rapidly consumes available 
memory as the size of the problem increases. Due to memory 
limitations, the algorithm Held and Karp presented is exact 
only up to 13 cities: They also present a dynamic progran- 
ming method which approximates the solution for larger prob- 
lems. The approximation is done by partitioning the large 
problem into several smaller problems or subproblems, each 
subproblem containing 13 or fewer nodes. Each partitioned 
subproblem is solved for its optimal tour resulting in a set 
of suktours for the original problem. The subtours are then 
broken and linked to each other to provide a tour for the 
original large problen. 

Little et al. (Ref. 6] present a branch and tound 
method. The algorithm branches on whether a particular arc 
is included or excluded from the tour being constructed. 


This kranching partitions the solution set containing all 


tours into subsets. Each arc in the total graph can be 
associated with some subset of tours, and the idea is to 
find the subset containing the optimal tour. For each 


subset of tours, a lower bound is computed foc the tours 
within the subset. The tour is constructed aS arcs with 
favorable lecwer bounds are included in the tour. The algo- 
rithm continues branching and computing lower bounds until a 
subset is found which contains one tour and this tour's 
distance is less than or equal to the lower bounds of the 
other sufksets of tours. This tour will be the optimal solu- 
tion to the TSP. The memory reguired for this technigue may 
re extensive and computing time grows exponentially with the 
size cf the problen. 

Many tour construction heuristics have been proposed. 


Rosenkrantz, Stearns, and Lewis [Ref. 2] and Golden et al. 
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[Ref. 4] review several of these heuristics. [The two most 
popular forms of these heuristics are the nearest neighbor 
algorithm and the insertion technigues. 

The nearest neighbor algorithm operates exactly as its 
name implies. A starting node is chosen and its nearest 
neighkoring node is selected as the next node in the tour 
keing constructed. The algorithm then iteratively selects a 
node not yet in the tour which is nearest to the previously 
selected node to enter the tour. This selection continues 
until all nodes have been selected and then the last node 
selected is connected to the first node to complete the 
tour. 

The insertion technigues include the nearest, farthest, 
cheapest, and arbitrary insertion algorithms. Insertion 
technigues begin with a starting node and select the nearest 
node to create a subtour. The algorithm iterates selecting 
a node which is nearest to or farthest from any other node 
in the subtour, or arbitrarily selects the next node to 
enter the tour. This selected node is then inserted into 
the subtcur wherever it provides the least increase (or in 
the case of the farthest, the greatest increase) in distance 
in the new subtour. These heuristics generally providea 
suboptimal solution to the TSP, but they have the advantage 
of being very simple to implement and have very fast solu- 


tion times with minimal memory consumption. 


Be SUBTOUR ELIMINATICN 


The Subtour Elimination solution method [Ref. 1] attacks 
the ISP by solving an nby naSssignment problem with the 
added constraints that the final solution must contain a 
cycle and the cycle cannot be of size n-1 or less. The 
number cf integer constraints necessary to eliminate all 


subtours or cycles of size n-1 or less iS enormous. 
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Therefore, the initial solution is generally found using a 
trelaxedt form of the aSSignment problen. The relaxed 
assignment problem cmits the subtour constraints in its 
formulation. 


Mathematically stated the relaxed assignment problem is 


qzt 37! 
Sut. - cae i= 12.2 
j- 
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where xi, 0708 X zi =1 


The solution to the relaxed assignment problem provides 
the initial lower bound on the optimal tour length. If the 
optimal assignment sclution is a tour then it is an optinal 
soluticn to the TSP. If the assignment solution is not a 
tour then there exist subtours which must be removed until a 
Single tour exists. If a subtour exists with k arcs in the 
subtour, then there are k possible subproblems to which the 
Eroblem may branch. These k subproblems each have an addi- 
tional constraine, the constraint for subprotlen i 
(l=1,...,k) being the exciusion of arc 1 from the problem to 
eliminate the subtour. Instead of branching into k subrrob- 
lems the aprroach is to branch into two subproblems. An arc 
is selected from the subtour and the exclusion of this arc 
becomes a new constraint in one of the subproblems while 
1nclusion of the arc in the tour becomes a nhew constraint 
for the alternate subproblen. The modified assignment 
problem is solved again for each of these new constraints 
and provides the lower bound for its respective branch. In 
this way only the necessary subtour elimination constraints 
are added as needed to the assignment problem, rather than 
attempting to add all possible subtour constraints to the 
initial assignment froblen. Ideally the branch solution 
subset with the lowest bound is then selected to be tested 


Apo 


moesee iieat is a tour, and, ifeiee 2s not, the prcblen 
Eranches again on the lowest bound. iepes Joba syovel | shine 
continues until all subtours have been eliminated and a tour 
is found which has a distance less than or equal to the 
lower bounds of the cther solution subsets. 

There are two kEasic methods for searching the tree 
created Ey the branch and bound process. One method is a 
breadth-first search. As the program is branching and 
computing bounds, a tree is created with leaves or terminal 
nodes. A breadth-first search reguires storage for the 
soluticns of each leaf and a search through these leaves for 
the next subproblem on which to branch. This requires 
extensive, usually exponential, storage. Branching at the 
lowest bound requires a breadth-first search. The second 
method for searching the tree is a depth-first search. The 
depth-first search is a more localized approach to handling 
the search. The depth-first search need only store the 
immediate sclution while searching a branch. If during the 
search it is found that no further branching is possible 
from the current node, the process must return to a previous 
leaf to continue the search for the optimal tour. The 
depth-first search kacktracks and recreates the previous 
solution at that leaf. The depth-first search does not 
require exponential storage like the breadth-first search, 
but it usually requires more computation time. 

The subtour elimination method is described as Leing 
exact. The subtour elimination method usually requires less 
computation time than the branch and bound approach of 
little et al. since the bounds obtained with the assignment 


relaxaticn are tighter. 
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C. TEE ONE-TREE FORSULATION 


The one-tree formulation was introduced and developed by 
Held and Karp [Ref. 7,8]. The general method for the solu- 
tion of the TSP begins with the construction of a mininun- 
weight ‘cne-tree.* if the one-tree 1S Not (a tour Stneaeay 
integer linear program may be used to oktain the final solu- 
elon 

Io assist in the description of this solution method 
several terms need defining. First a tree is defined as a 
connected graph without cycles. A minimum-weight spanning 
tree is a tree with finimum total weight on its edges. A 
one-tree is a tree to which one edge has been added yielding 
exactly one cycle. A minimum-weight one-tree is a one-tree 
with minimum total weight on its edges. Finding a minimun- 
weight one-tree is a relaxed version of a TSP since the 
soluticn to a TSP is a minimum-weight one-tree aAaving every 
vertex of degree two. 

In the solution method proposed by Heli and Karp 
[Ref. 7] the algorithm begins by constructing a minimun- 
weight one-tree. A Mminimum-weight one-tree can be can be 
found by first constructing a minimum-weight spanning tree 
through noces 2 to nand then adding to the graph the two 
arcs of least weight from node 1 {[Ref. 9]. A minimun- 
weight one-tree may also be constructed by constructing a 
Minimum-wéeight spanning tree through nodes 1 to no and then 
add to the graph the next minimum-weight edge not yet used 
from the distance matrix. If the one-tree constructed is a 
tour (i.e. the one-tree's vertices are all of degree two) 
then the TSP is solved. Otherwise, the one-tree must be 
converted into a tour. 

Integer linear programming is used to obtain an optimal 
tour if the one-tree formulation does not satisfy a tour. 


Held and Karp introduce the concept of a ‘gap* function and 


21 


then use a special integer linear projram to minimize this 
function. In essence there is a transformation on the one- 
tree variables which allows the integer linear program to 
search for the optimal tour. Held and Karp suggested two 
methods for optimizing the integer linear program, an ascent 
method and a branch and bound method with the ascent method 
embedded in it. Other improvements to the Held and Karp 
method were suggested by Hansen, Krarup [Ref. 10] and Houck, 
Picard, Vemuganti [Ref. 11], Bazarra and Goode [{Ref. 12] 
did further work on Held and Karp's approach usiog optinmiza- 
tion of a lagrangian dual in lieu of solving integer linear 
program. They proposed a branch and bound scheme with subg- 
radient optimization of the dual to transform the one-tree 
into the optimal tour. The one-tree method as described is 


exact and may reguire extensive computation time. 


De. TCUR IMPROVEMENT 


The Tour Improvement method assumes that an arbitrary 
tour is available. The method operates by perturbing or 
exchanging arcs in the tour until a better tour is found. 
The Tour Improvement method terminates when a better tour 
cannot be found. In its simplest form, it is possible for 
this method to check (n-1)!/2 arc exchanges for the symnme- 
tric ISP before terminating (i.e., it investigates all 
availatle answers). 

Dantzig, Fulkerson, and Johnson [Ref. 13,14] proposed a 
algorithm which starts with an arbitrary tour and uses 
integer linear programming to improve the tour. Theslse Vs 
transformed into an integer linear program and solved using 
the simplex method. DantZig, Fulkerson, and Jchnson 
[Ref. 15] mention that the number of constraints needed to 
characterize the problem is ‘astronomical. Instead of using 


all the constraints, they begin by using a relaxed version 
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of the problem and add constraints as needed to™the integer 
linear program to maintain feasibility. In their method 
infeasibility appears as ‘loop* oor subtour. To continue 
solving the problem a constraint is added which removes the 
subtour yet does not eliminate any of the available tours. 
This algcrithm reguires the addition of constraints during 
the soluticn process to maintain feasibility. Several 
researchers have suggested improvements on Dantzig, 
Fulkerson, and Johnson's method for solution of the TSP. 
(Ref. 16,17,18,19] 

Croes [Ref. 20] proposed a heuristic whereby a subset of 
all possible transformations is tested, a transformation 
keing the transformaticn of one tour into another tour. He 
called these transformations ‘inversions because they 
inverted the sequence of nodes in part of the tour to create 
a possible improvement. Lin [Ref. 21] later describes Croes 
‘inversicn free tours' as 2-optimal tours and goes on to 
describe k~optimality where k is some fixed number less than 
ne The general idea is to exchange k arcs iteratively in 
the tour while testing for improvements. Lin and Kernighan 
[Ref. 22] presents a modified k-opt algorithm where k is not 
fixed. In this algorithm k May vary Eronm 2 ten ie The 
floating k-opt algorithm is also described by Christofides 
and Elon her soar! 

The advantage of the ‘inversions' or k-opt solution 
method is that for a given problem, the memory needed to 
solve the problem is fixed. There are no constraints which 
must ke added to the froblem, and the decision rule used is 
Simple. The k-opt method is exact, but may be terminated 
before optimality is reached to provide a satisfactory 
(suboftimal) solution. This generally reduces run time and 
still provides a gocd solution to the problem. The k-opt 
method may also be programmed to reyuire very little user 


interaction; in the best case the user need only enter the 


2S 


distance matrix for the nodes desired in the tour. The 
k-opt algorithm may require extensive computation time, 
particularly if during solution k is found to ke greater 
than 2 on many of the iterations. 

This chapter has presented a general survey of some the 
solution methods for the TSP and some algorithms associated 
with these solution methods. The next chapter will discuss 
which of these solution methods were employed to solve the 
bouytender problen. The buoytender routing algorithm will 


then Ere presented and discussed in detail. 
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IV. BUOYTENDER ALGORITHM 

This chapter will discuss the algorithm with which we 
propose to solve the buoytender problen. The reasons fen 
the algorithm's selection and some background on its devel- 
opment will be presented, together with a detailed descrip- 


tion of how the algorithm addresses the buoytender problen. 


Aw. ALGORITHM SELECTICN 


The algorithm is a combination of two of the previously 
mentioned solution methods in that it combines a fTour 
Building heuristic with a LOuGe ate Tour Improvement 
heuristic. The ccncept of combining a Tour Building 
heuristic with a 2-cft or 3-opt heuristic was proposed by 
Golden et al. (Ref. 24] as a relatively fast computational 
solution method which would provide an optimal or near 
optimal tour. Their ‘composite! alyorithm is the foundation 
on which the buoytender algorithm was developed. 

The heuristics used in the buoytender algorithm are the 
nearest neighbor algorithm and the k-opt algorithn. The 
nearest neighbor algorithm is used to construct an initial 
tour while a versicn of the k-opt proposed by Lin and 
Kernighan [Ref. 22} is used to improve this initial tour. 
The nearest neighbor and k-opt heuristics were selected for 
use in the buoytender algorithm because very little interac- 
tive work is required, the interactive work is simple, no 
constraints need to ke added to maintain feasibility, comfu- 
tational time is reasonable, and the algorithm operates with 
a fixed amount of temory. 

The kEuoytender algorithm as presented requires very 


little user interaction. The nearest neighbor algorithm 
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needs only the entry of the distance matrix for the prcblen 
and the selection of a starting node. The k=-opt algorithn 
needs the same distance matrix andan initial tour. The 
heuristics! decision processes will provide a feasible and 
possibly optimal solution ina reasonable period of time. 
The nearest neighbor algorithm was selected to construct the 
initial tour because it iS computationally faster than 
insertion technigues, and provides a good initial solution 
to the problem. [Ref. 25] 

A critical factor in solving the ISP on a microcomputer 
is the memory reguired by the algorithm used. The amount of 


memory reguired for solution of a TSP of size n is fixed for 


the nearest neighbor and k-opt heuristics. With known 
memory reguirements for the heuristics, a program can be 
developed tc fit the Coast Guard Standard Terminal. The 


nearest neighbor algcrithm reguires an n by n distance 
matrix and an nm arrary in which to store the tour as it is 
constructed. The k-cpt algorithm requires the same storage 
as the nearest neighkor plus an by n decision matrix, anda 
2n array for recording nodes selected for inversion in the 
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Be. PROGRAMMING THE ALGORITHM 


The Euoytender algorithm was developed as a_ series of 
modules which were then programmed as subroutines in the 
ecperating program. The algorithm was programmed in FORTRAN 
for operation on the Coast Guard Standard Terminal. FORTRAN 
was selected because it is efficient and a majority of 
programmers are familiar with it. 

There are some differences between the heuristics in the 
literature and the buoytender algorithm. The nearest 
neighkor algorithm as used by the buoytender algorithm is 


the same as outlined by Rosenkrantz, Stearns, and Lewis 
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fRef. 2}. The k-opt heuristic used in the buoytender algo- 
rithm is a modified version of the heuristic presented by 
Lin and Kernighan [Ref. 22]. There are three major differ- 
ences between the k-cft in the bouytender algorithm and the 
k-opt suggested by Lin and Kernighan. First, the Lin and 
Kernighan k-opt heuristic selects all possible nodes which 
will result in a sherter tour before creating the improved 
EOUT s The buoytender algorithm improves the tour by 
exchanging nodes whenever a favorable selection is found. 
This improvement-as-ycu-go procedure results in a simplifed 
selection process for the next set of nodes to be tested. 
Second, the Lin and Kernighan k-opt heuristic has an added 
facility for ‘limited backtracking' used when a particular 
exchange gives an improvement of zero. The backtracking 
procedure searches until a gain greater than zero can be 
found which imprcves the tour. In the interest of short- 
ening computational time, the buoytender algorithm does not 
‘backtrack! when a gain of zero is found, but instead treats 
zero gain aS no improvement and continues with the next 
selection. The third major difference between the Lin and 
Kernighan k-opt heuristic and the buoytender algorithm is 
related to what Lin and Kernighan call ‘reduction.* Ifews 
heuristic, after producing several locally optimal tours, 
checks fcr arcs which appear in each of the tours. These 
‘good' arcs are not allowed to be broken in further computa- 
tions for other locally optimal tours thereby reducing the 
humber of links to be checked for improving the tour. While 
decreasing run time, this procedure reguires more memory and 
is therefcre omitted from the buoytender algorithn. 

The ruoytender algorithm can be divided into three kLasic 
operations. First, the algorithm selects a random sample of 
k nodes, k<n. This random sample of nodes becomes the set 
of starting nodes for the next major operation in the algo- 


rithn, to create k nearest neighbor tours. The third 
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operation is to use the k nearest neighbor tours as initial 
tours for the k-opt improvement. The best of the k-opt 
improvement tours is presented as the solution to the 
problen. 

A further modification was made to the algorithm based 
on some computational results obtained from use of an early 
version of the algorithn. During programming of the kuoy- 
tender algorithm several benchmark statistics were cbtained 
from the pregran. Two of these statistics were the initial 
tour distances produced in the nearest neighbor phase, and 
the distances of the final tour solutions produced by the 
k-opt phase of the algorithm. There seemed to bea rela- 
tionship between the tour froduced by the nearest neighbor 
phase and the improved tour produce by the k-opt phase. A 
least-squares linear fit of the initail tour distance vs. 
the final tour distance was done for four data _ sets. In 
each case the fit prcduced a positive slope. This implies 
for a relatively large nearest neighbor tour, the final tour 
from the k-opt phase will be relatively large. A relatively 
small nearest neighbcr tour will produce a relatively small 
tour frem the k-opt phase. Using this information the Luoy- 
tender algorithm was further modified to reduce computation 
time without compromising too much of its ability tc produce 
Optimal] solutions. For a set of x nodes, k<n, the nearest 
neighbor algorithm provides k initial tours. The set of k 
nodes is selected at random from the population of n nodes. 
Then for a set of L nedes, L<k, the k-opt algorithm is run 
Homo@oucain L £inal tours. The set of L nodes are drawn fron 
the k nodes which froduced the L shortest initial tours. 
The shortest final tour from the L final tours produced by 
the k-opt phase will be the algorithm's solution to the 
problem. The final version of the algorithm has a greatly 


reduced operating time over the initial version. 
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C. ALGORITHM DESCRIPTION 


In this section, a detailed description of the buoy- 
tender algorithm is given. In this description references 
will be given to the program subroutines so that the reader 
May associate the portion of the algorithm under discussion 
with its operational counterpart in the progran. 

The kuoytender algorithm begins by reading a file from 
mass storage containing the number of aids (nodes) in the 
problem, the distance matrix, and the names of the aids 
(subrcutine OBTAIN). A program MATUTL.FOR, which assists 
the user in creation of this data, can be found after the 
listing ef the kuoytender route selection progran, 
ROUTE.FOH. 

After oktaining the data the algorithm selects a randon 
sample of nodes, baséd on the size of n, from which it will 
create tours (subroutine NODE). The function which selects 


the size cf the randcm sample is 
Sample Size = 5 + [2( Log{ n-4 ))\. 


The function is designed to capture almost the entire popu- 
lation for testing when the problem 1S small and to capture 
only a small portion of the population being tested when the 
problem is large. This produces a sample size which 
increases the probakility of obtaining the glokal optimun 
when the problem is small yet produces a sample size which 
will maintain reasonakle computation times when the problem 
is large. 

After determining the sample size, a random sample of 
nodes is taken. Since there is no random number generator 
for FORTRAN installed on the Coast Guard Standard Terminal, 
a pseudo-random number generator was added to program (func- 
tion RANDOM). The coding tor the generator comes from 
Wichmann and Hill [Ref. 26] ; it 1s described by “abe 


oe 


authors as an efficient psuedo-random number generator 
having a cycle length of 2.78 x 1013 . Nodes are sampled at 
random without replacement from the integer population of 1 
to n for use as starting nodes for the nearest neightor 
phase. 

The algcrithm then moves into the nearest neighbor fhase 
(Subroutine NBR). A tour is created for each node in the 
randcm sample, with the node from the sample as the starting 
node of the tour. The algorithm selects the nearest node to 
the starting node for the next entry in the tour construc- 
ION « The algorithm then iterates selecting a node not yet 
in the tour which 1s nearest to the previously selected node 
as the next entry in the tour. This process continues until 
all nodes have been selected, at which time the last node is 
connected to the starting node to complete the tour, and the 
tour distance is computed. After all the initial tours have 
been constructed they are ordered by their total distances 
from shortest to Ilcngest. Tie wom SHOTLCS Ie TOUrS are 
retained for possible improvement in the k-opt phase. 

A deécision matrix is created which designates waich 
distance arcs in the distance matrix are presently members 
of an initial tour (subroutine MARKD). The algorithm then 
proceeds to step through each node of the tour testing it 
moe kK-Optimality. 

The tour 1s ‘prepared* at each node to be checked for 
k-optimality (subroutine TRPREP). This process simplifies 
some of the later operations. Hetejwmeekt. Iseaheiprulea to 
consider the tour as a circle with n postions on its circun- 
ference (see Figure 4.1), and to let positions on the 
Circular tour be denoted by p(i) where 1=1,...,n. A node in 
the tour is represented as a position on the circle. MTRPREP 
places the starting node in p(1) and places the remaining 
nodes in the tour at p(2) through p(n) where (1) 1s the 


node's pesition in the tour. 
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Figure 4.1 Circular Representation of a Tour 


Since floating k-cpt algorithm is in essence a series of 
2-opt iterations a description of the simpler 2-opt prin- 
ciple will ke presented first. Figure 4.2 shows the 2-opt 
principle in graphic detail. In Figure 4.2 (A), the tour is 


represented in circular forn. The tour 1S now perturbed to 
try and find an improvement. To perturb the tour inte 
2-opt procedure, twe arcs will be broken in the existing 


tour and two new arcs will be selected as incoming to form a 
new tcur. This provides a gain function which can be used 
to check for improvement. The gain function is 

G = E (broken arcs) - = (incoming arcs). 
If G ig fositive this set of arc exchanges will improve the 
tour. 

Tc perturb the tour, one of the arcs incident to node 
being checked for 2-optimality must be broken or removed 
from the tour. In the buoytender algorithm the node feing 
checked for k-optimality will always be in p(1) of the 
circular tour andthe arc being broken will always be 
retween p(1) and p(2) of the circular tour. An incoming arc 
must then re selected to replace the broxen arc. This 
incoming arc must ke selected from those arcs with a 


starting node which is the same node as the end node of the 
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Figure 4.2 Graphic Description of 2-opt 


previously kroken arc. Figure 4.2 (B) Shows scme of the 


available selections for the incoming arc. The incoming arc 


BZ 


is in general the shortest available arc, but may be any of 
the available arcs which maximize G. This subtour must be 
troken and another incoming arc must be selected so as to 
reconstitute a tour. Therefore, one of the tour's arcs 
incident to the end node of the first incoming arc must be 
troken. The end node in Figure 4.2 (C) is at p(7). There 
is only cne of two possible tour arcs which may be broken to 
allow the incoming arc to become part of the tour. Figure 
4.2 (D) shows the selection of the wrong incident tour arc 
for removal. The creaking of this arc and the subsecuent 
reconnection to the ncde being checked for 2-optimality does 
not result in a tour. Instead, two subtours are created 
giving an infeasible situation. Figure 4.2 (E) shows the 
Freaking of the proper incident arc and subsequent reconnec- 
tion te the original node to create the ‘improved tour’. 
Also notice that some of the arc directions in the original 
tour must be reversed to provide a consistent direction of 
movement through the tour. Reversing these arc directions 
is the same as inverting a seyuence of nodes in the tour. 
The result is shown in Figure 4.2 (F) as anew ‘timproved' 
CLEC UlaAE Coun. 

The 2-oft procedure is carried out using each node in 
the tcur as the starting node checked cor 2-optimality. One 
interation of the 2-opt procedure through all n nodes will 
produce a tour which is 2-optimal. 

The K-opt algorithm is very similar to the 2-opt proce- 
dure descriked above. Instead of just reconnecting the tour 
as the 2-opt procedure does (p(6) to p(1)), the k-opt checks 
to see if another incoming arc, starting at the end nede of 
the last broken arc, can be selected to further improve the 
tour (Figure 4.3). The procedure interates until no further 
improvement can be fcund. A k-opt, using an example with k 
equal to four, might look like Figure 4.4 The floating k-opt 
algorithm has the ability to vary k during the k-opt 
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Figure 4.3 K-opt Selecting Next Incoming Arc 
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Figure 4.4 Example of 4-opt 


solution process. This is more desirable than fixing k in 
the solution process. For example, if k were fixed at four, 
the algorithm can only check four are exchanges for improve- 
ment but the problem may need a five arc exchange to produce 
the optimal solution. By allowing the k-oft to vary k, the 
probability of the algorithm finding the optimal solution is 
increased. The floating K-oOpe is bounded between 
2-optimality and n-optimality. The k-opt cannot be greater 
than n-optimal because an n-opt exchange creates the optinal 
[SOUL . 
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The k-opt procedure in the buoytender algorithm operates 
by checking each node in a tour for k-optimality. The algo- 
rithm, using subroutine TRPREP, creates the circular tour 
for the node being checked for k-optimality. The circu@ac 
tour is prepared so that the initial arc being broken will 
always be between p(1) and p(2). The algorithm checks the 
circular tour at the starting node for k-optimality. If no 
improvement can be made breaking the first are between p(1) 
and p (2) then the circular tour's direction is reversed 
(subroutine REVTOR) so that the alternate incident arc is 
placed between p(1) and p(2) of the circular tour and this 
tour is checked for k-optimality. 

The kreaking of the arc between p(1) and p(2) requires 
an inccming arc to replace it (subroutine SELCTY). The 
algorithm selects the five shortest available arcs which 
originate at position two of the circular tour and terminate 
at ancther position on the circular tour. The arc from p (2) 
to p(s) is not available since it is already in the tour, 
and the arc from p(2) to p(1) 1s not available since it is 
the arc just broken to improve the tour. For each of the 
selected incoming arcs another arc must be broken to elimi- 
nate a subtour (as in the 2-opt procedure). It is then 
necessary to reconnect the arcs to reconstitute the tour. 
This results, for each of the selected arcs, a set of two 
1ucoming arcs and two outgoing arcs. Using the previously 
defined G, each of the five sets of arcs is checked to see 
which maximizes G. If G>0 then the exchange or inversion is 
made (surroutine ADJTCR), otherwise no exchange is made. 

If an exchange is made the algorithm goes back and 
searches for a new incoming arc from the end point of the 
last Eroken arc. Sukroutine ADJTOR has inverted the tour so 
that the last arc used to reconnect the tour is now between 
p(1) and p (2) of the circular tour. Subroutine SELCTY is 


again used to select the next best available incoming arc as 
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described atove. The algorithm interates through selecting 
an incoming arc and inverting the tour as long as improve- 
ment, 1.@€. foSitive gain, can be found. 

When improvement can no longer be found the direction of 
the tour 1s reversed at the last node where improvement was 
found, and the tour is tested again for possible further 
improvement. When nc further improvement can be found and 
the reversed tour has also been tested, the algorithm then 
increments to the next node in the tour and starts the k-opt 
procedure again testing this node for k-optimality. Each 
node in the tour is used as a starting node in the k-opt 
procedure. When all nodes in the tour have been tested for 
k-optimality the k-opt phase ends and the resulting tour is 
stored. 

After each of the five shortest initial tours have been 
tested for k-optimality, the shortest of the improved k-opt 
tours 1s selected as the solution to the buoytender frotlen. 
Figure 4.5 is a flowchart of the buoytender algorithn. 

The next chapter discusses the results of some test 
problems solved by the algorithm and makes’ concluding 
remarks about the algorithm. 
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V. PROGEAM RESULTS AND CONCLUSIONS 
This chapter will discuss some results obtained from 
uSing the kuoytender route selection progran, including 
results from an actual problem from the U.S. Coast Guard 
Cutter BLACKHAW. All program results were obtained using 
the Coast Guard C3 Standard Terminal. The chapter will 
conclude with some observations regarding the bucytender 


procklem, the algorithm, and the test results. 


Aw. FEROGKAM RESOLTS 


To program and test the buoytender algorithm it was 
necessary to obtain problems with known optimal solutions. 
Initially, avery simple and small problem with nine nodes 
was Created. This test problem's optimum tour is the most 
obvious tour obtained by inspection of the problem since the 
nodes are arranged in an almost circular pattern. 
Connecting these nodes following this nearly circular path 
results in the optimum solution. Several alternate tours 
were checked by hand computation and were found to be subop- 
tinal. Therefore, the near circular path for this test 
problem is assumed to be optimal. When the program was 
tested with this test problen, the program produced this 
hypothesised optimum solution with a computation time of 13 
seconds. 

With the knowledge that the program appeared to be oper- 
ating correctly on the initial test problem, further testing 
was necessary to check the program's ability to solve other 
problems. Since the Euoytender route selection problem is a 
symmetric case of the TSP, the literature was searched for 
test protlems. Three articles [Ref. 5,20,27] yielded eight 
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symmetric TSP's of various sizes. These eight literature 
Eroblems ranged in size from 5 to 57 nodes with known 
optimal answers. The program, when tested with these frob- 
lems, produced solutions which ranged from the optimal to 
within 3.5 percent of the oftimal distance. 

Another problem was created by the random selection of 
70 pcecints on a Euclidean plane. A distance matrix was 
computed for these 70 nodes and this problem was solved by 
the program. The problem's optimal solution is not known, 
but the problem was run to ensure the program could handle 
the maximum problem for which it was designed, and to 
produce a sample solution time for this maximum sized 
problem. 

With evidence that the program was solving test prob- 
lems, a further test was conducted using an actual kuoy- 
tender routing problem, with data provided by the U.S. Coast 


Guard Cutter BLACKHAW (WLB-390), homeported in San 
Fransiscec. The aids to naviyation run data was for a trip 
scheduled from 9 April 1984 to 16 April 1984. With the 


experience and historical data this vessel has availatle 
from servicing these aids for many years, it is hypothesized 
that the route which was scheduled for this ATON run would 
be optimal or nearly optimal. The planned route for the 
ATON run had a distance of 455 nautical miles. When this 
problem was entered and solved by the buoytender program, 
the sclution was also 455 nautical miles. Although the 
distances were the same, the scheduled route and the progran 
route differed slightly in their sequence for visiting the 
aids. This difference may be attributed to the fact that, 
as stated in Chapter II, the triangle inequality may be 
reguired to be a strict equality on some of the arcs. This 
difference is possible since the buoytender may, due to the 
strict equality on some arcs, Le reyuired to ‘pass' a previ- 


ously serviced aid tc complete the route. Figure 5.1 shows 
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a five node example of this situation. This suggests there 
may be several sequences of nodes which will provide an 
optimal solution to the buoytender problen. 


TABLE 1 
Tours for BLACKHAW Problen 


Route as Scheduled Route Produced 
by the Buoytender by tne Program 
1 YB ISLAND 1 YB ISLAND 
Z MONTARA LWB10A 20 MOSS LBB MLA 
3 N ISLAND LWE8 19 MONTRY BY LB B 
4 SANTA CRUZ ME 18 MONTRY HBR MB 
5 PLANCAS LWB4A 17 PT PINOS LWB2 
6 SAN SIMEON LEB! 16 PT CYPRESS LGB6 
a, VON HELM R LGEY 5 P BLANCAS LWBYA 
§ ESTERO BY GB10E 6 SAN SIMEON LBB1 
9 SOUZA R LGB , VON HELM R LGB& 
10 WESTDAHL R LEBF!1 15 ESTERO BY LWB 
11 SANLUIS LWE3 8 ESTERO GBIOE 
12 LANSING RB 14 4ORRO BY LBB! 
13 ET BUCHON LWE2 13. PT BUCHON LWB2 
14 MORRO BY LBB1 10 WESTDAHL R LBB! 
15 ESTERO BY LWE 9 SOUZA R LGB 
16 ET CYPRESS LGE6 12 LANSING R B 
17 ET EFINOS LWB2 11 PT SANLUIS R B 
18 MONTRY HBR ME 4 SANTA CRUZ MB 
19 MONTRY BY LB E é' AN ISLAND LWB8 
20 MOSS LBB MLA 2 MONTARA LWBIO0A 
1 YB ISLAND 1 YB ISLAND 





Figure 5.1 Multiple Optimal Solutions 
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TABLE 2 
Results of Test Problems 


Problem & Sclucion TOUE DISt. (eves) Solution 
(Ref.) Size Time (Sec) Program Optimal Quality* 
27 5 5 148 148 100% 
test auth.) 9 ie 22205 Dos Ie 100% 
Barachet(2/) 109 17 376 378 100% 
3LACKHAA 20 50 455 455 * 
Croes (20) 20 52 254 246 10 3. 25% 
HeldKarge (5) 2 Ue 1711 1711 100% 
Kar (any See 121 109 29 10861 1002 e225 
Dan era 4 2 188 705 699 100. 863, 
HeldKarf(5 48 Zoo 11847 11470 103. 28% 
Kar (27f 57 340 133 6y 12956 103.17% 
Max (auth.) 70 494 6791.07 “unk: - unk 
* 100% is optimal; 103.28% is 3.28% greater than optimal. 
** The distance 1S hypothesized to be optimal. 
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Figure 5.2 Least Squares Fit of Size vs. Solution Time 


Table 2 presents a summary or the eleven test precrlems 
Showing proklem sizes, solution tines, and the guality of 
the solutions produced by the program. All of the run times 
are reasonatle by the criteria stated in Chapter II. The 
largest prcblem, 70 nodes, had a solution time of 494 


seconds, or roughly a little more than eight minutes. 


43 


Ccmputation times for the eleven problems were plotted 
against froklem size, as shown in Figure 5.2 A least-squares 


prolynemial fit of the data points is 
Time = .584 + .62n + .0914n2 + .00000994n3 


and this curve is alsc shown in the figure. 


B. CCNCLUSIONS 


The kbuoytender route selection algorithm provides 
another useful way of scheduling ATON routes. As a tool for 
the operator to assist in planning the route to be taken, 
the program should produce optimal or near optimal solutions 
for precbtlems up to a Size of 70 aids. The user may accept 
this route aS 1s or may desire to modify the route kased on 
circumstances requiring human judgement. Appendix A isa 
user's guide for the matrix utility and buoytender route 
selection programs. 

The Euoytender route selection algorithno, with the 
nearest neighbor and k-opt heuristics, provides a guick and 
Satisfactory soluticn to a symmetric traveling salesman 
problem cn problems up to a size of 70 nodes. The problen 
may have several rceutes with different sequences for 
visiting the aids, tut the routes may all have the sane 
distance. 

It is hoped that this algorithm will be of use to the 
fleet cf U.S. Coast Guard buoytenders. As a decision aid it 
May help the operators of buoytenders obtain optimal or 
hearly optimal routes to service their aids, possibly saving 
time and fuel. 
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APPENDIX A 
USER'S GUIDE TO PROGRAMS 


The intention of this appendix is to provide instruc- 
tions for the use of the buoytender route selection program 
and utility programs. It should be noted that pricr to 
using the Luoytender route selection program, the user will 
need to use the matrix utilities program listed at the end 
of this appendix to frepare a data file. 

The kuoytender route selection program and the matrix 
utilities progran are written in Convergent FORTRAN. 
Convergent FORTRAN is FORTRAN 77 compatible. The progrags 
should be typed in as presented and linked to the operating 
system. The Convergent FORTRAN manual has instructions 
concerning compiling and linking. Recommended names for the 
run files are ROUTE.RUN and MATUTL.RUN. Both programs will 
operate cn either a IWS or AWS station with 264K memory or 
more. The programs run about three time slower on an AWS 
terminal than on an IWS terminal. 

The user will be reguired to enter text and numerics 
into the programs. All numeric entries are to te integer 
entries (no decimal point) with the exception of the 
distance entries. All distance entries require a decinal 
foint entry for proper input. 

MATUTL. RUN is capable of entering problems up to a size 
of 100 aids. ROUTE.RUN is designed to handle problems up to 
a size of 70 aids. MATUTL.RUN allows the entry of a large 
matrix of aids and includes a utility which allows the user 


to select a subset of this large matrix for use in the 


program ROUTE.RUN. This allows the user to type in one 
large matrix and tken create smaller matrices from the 
larger as needed. This feature will preclude the user fron 
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typing in a distance matrix every time he desires to run a 


problen. 


A. USE CF MATUTL. RUN 


The matrix utilities program is designed for creating 
and editing data files which will be used in the buoytender 
route selection progran. A data file will contain the size 
of the problem n, the distance matrix for the problem, and 
the names of the aids to be visited. 

The size of the froblem, n, will be the number of aids 
to be visited plus one. The plus one is for the port fron 
which the kEuoytender will start and end the ATON run. 

To create the distance matrix some extensive plotting is 
reguired. Most buoytenders have a set of tracklines they 
regularly follow to visit their aids. In nost cases these 
existing tracklines will be sufficient to use in the 
distance matrix. To create a distance matrix for entry into 
a data file it is reccmmended the user use the format shown 
in Figure A.1 . Distances should be computed in nautical 
miles and may range from .01 nm to 9999.99 nn. Figure A.1 
is an example of a prcblem of size nine. 

Since the problem 1S symmetric, the user need only 
calculate and enter the distances as shown in Figure A.1. 
The matrix utilities will automatically complete the rest of 
the matrix entries. The distance you should compute should 
be the shortest possible distance from the row entry to the 
column entry. Using the example in Figure A.1 the entry in 
row 3, column 6 is the shortest possible trackline frem aid 
neo aid £E. Once tke entries for the distance matrix have 
been prepared the user is ready to use MATUTL.RUN. 

MATUTIL.RUN has ¢cme Simple menus and some descriptive 
comments to helfpr the user. The menus reguire a numeric 
response for the selection of an option. The first menu 


presents the utilities available to the user, Figure A.2. 
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Figure A.1 Sample Distance Matrix 


_—— ee 
1 = create 4 distance na. eux 
2- extract a small matrix from a larger matrix 

3 =  Correet sa VEU Sl h Stored matrix. 
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6 - exit this prcgraa 

What is your selection? | 


Figure A.2 MATOTL.RUN Main Menu 


1. Create a Distance Matrix 
This utility is for the initial entry of a distance 
Matrix into wandata ta. 1ce The user will be asked 


Size cf the problem, the homeport and aids names, 
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for the 
and the 


distances. The size of the problem is the number of aids 
plus one, for homeport. The utility will next ask for names 
of homeport and the aids to be visited. These names will be 


associated with a row and column number of the distance 


matrix. The nameS are text entries up to a maximum of 15 
enatacters. 
The utility will next ask for distances. It will 


ask fer the distance from a row number to a column number. 
If the user sets up the distance matrix as shown in Figure 


A.1 he will be able to then just read each row across its 


columns while entering the data. These distances must be 
real numckers. The distances should be nautical miles and 
May cange fron .017 nm to 9999.99 nn. Since these entries 


are real numbers there must be a decimal point in each 


entry. If, for example, the user needs to enter two 
hauticai niles, he should enter it as 2. <return>. Be 
careful with these entries; if a decimal point is not 


entered the distance will be off by a factor of 100. FOr 
guickly and easily entering data the numeric keypad on the 
keyboard may be used and <next> may be used in flace of 
<return> to enter the data. If an error is made while 
typing in numbers and <return> or <next> have not yet been 
pressed, then the error may be corrected by using the kack- 
space key and retyping the number. If an error is made and 
<return> or <next> have already been pressed, then the user 
Should ncte the row and column number of the error and use 
the correction utility to correct the entry. 

When the user finishes typing all entries the 
program will prompt fer a filename in which the data will be 
stored. This name may be up to 15 characters long. Jie opm 
some reason the user needs to end the terminal session 
before the entire matrix is entered, data entry may be 
stopped by typing ina negative distance. The user will 


then be asked for a filename for the partial data file. To 
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complete the partial data file the user must use the correc- 


iLO n (Utila tye 
2. Extract a Matrix 


This utility is designed so tnat the user may choose 
a subset of aids frem a larger matrix. The utility will 
create the proper data file for a subset of aids. 

The user should have an idea of all possible aids 
which may te visited on a general buoytender run (i.e. the 
spring run, the north run, the inside cun, etc.). The user 
should create a distance matrix using the create utility of 
all aids which could possibly be visited on a general run. 
Usually cnly some portion of these aids will be visited on 
an actual run (for example on a particular spring run only 
90% of the aids which might be visited on this run need to 
ke visited). In this case the user creates a matrix cf all 
aids which might be visited on a run and then each time this 
run is scheduled the user can then use the extract utility 
to guickly and easily create a data file for those aids 
which need to be visited. 
user first use the list utility on the larger or initial 
data file. This will provide the user a list of aid nanes 
in the data file and present them with their corresponding 
row number. 

The extract utility will prompt the user for the 
hame of the data file containing the general run'ts aids. 
The user is then asked for the size of the new eroblen 
wherelky the user enters some m, m less than n in the initial 
data file. The utility then asks the user m times for row 
humbers from the initial matrix. After m row numbers nave 
been entered the utility creates the new data file and askes 


the user for a filename for this new data file. 
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SemecORLect sa Previousiy Stored, Natrix 


This utility is to correct entries ina distance 
Matrix or tc complete a partially entered data file. The 
utility has a menu which askes if the user wishes tec make 
Single entry or seguential entry corrections. The single 
entry mode asks the user for the row number and column 
number of the entry desired to be corrected. The user then 
is asked for the new distance and is returned to the correc- 
tion menu where he may select to make another correcticn or 
exit the utility. 

The sequential corrections are Similar to the entry 
of distances in the create utility. The user is asked for 
the starting row number and starting column number where 
correcticns are to begin. The utility asks for the 
remaining column entries on the starting row and then will 
proceed to the next row until entries are terminated. The 
utility automatically terminates at the end of the wtatrix or 
the user may terminate the Sequential mode early by entering 
a hegative number. If a negative number is entered it will 
not appear in the distance matix. Terminating the sequen- 
tial mode returns the user to the corrections menu for 
another correction or to exit the correction utility. 

Upon exiting the correction utility the user is 
asked for a filename for the corrected file. This filename 


may be up to 15 characters long. 


4. Frepare a Distance Matrix for Hardcopy 


This utility is for preparing a hardcopy of a 


distance matrix from a stored data file. The data files are 


saved as seguential files, (iee., the data files area 
string of numbers). This utility will construct a file 
which will resemble a matrix. Due to the width and length 


limitaticns of the printer and paper, the matrix is printed 
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in sections which will have to be cut and pasted together to 
display the whole matrix. A hard copy of the whole matrix 
will provide the user with a means of checking to see if 
data entries are correct. If an entry error is found, the 
matrix also helps the user to locate the row and column 
humber of the error. 

The utility asks the user for the name of the data 
file from which a distance matrix is needed. The utility 
creates the hardcopy file andthen prompts the user to 
obtain a hardcopy by using the FORMAT command to print out 
file MATFRI.DATA. 


5- List Aid Names 


This utility provides a hardcopy of the aid names 
from a data file. This may be useful to the user when using 
the extract utility. The user is asked for the name of the 
data file from which the list shall be made. The utility 
then prepares a file containing the row numbers’ and their 
associated aid names. The user 1S prompted to use the 
FORMAT Cctimand to optaln a pEIntOUGeoreel Io. 01th. 

These matrix utilities are for assisting the user in 
creating and manipulating data files and distance matrices. 
A correc: data file is necessary for the proper operation of 


the buoytender route selection program, ROUTE.RUN. 


Be USE CF ROUTE.RUN 


The user of ROUTE.RUN is very simple once a_ data file 
has been created. The user runs the ROUTE.RUN program and 
he 1s prempted for the name of a data file. The user is 
then asked for the name of an output file for the results. 
The program obtains the size of the problem from the data 
file and will ccmpute an approximate solution time. This 


approximate solution time is based on results from some test 


om 


problems solved by the program and it assumes the user has 
version 8.01 of the operating system. If another version of 
the operating system is used the solution time aay differ. 
After displaying the approximate solution time the 
program reads the rest of the data file and proceeds to do 
the necessary computations. The user, depending on the size 
of the problem, will have just enough time to get a cup of 
coffee before obtaining the solution. The solution will be 
displayed on the screen and will also be placed in the file 
previously named for results, from which a hardcopy may be 
made. The user now has a possible route to take for that 


particular aids to navigation run. 


DZ 


A 
SAMPLE OUTPUT OF EUOYT 


PPENDIZ B 

ENDER ROUTE SELECTION PROGRAM 
The solution for data file BLACKHAW.MAT 
Total tcur distance = 455.00 


The tour is as follows 


1 YB ISLAND 

20 MOSS LBESaEA 

19 MONTDRY EY Tees 

18 MONTEY HBR MB 

17 PT PINOS LWB2 

16 PL CYPRESS EGes 

S P BLANCAS LWB4A 

6 SAN SIMEON LBB1 
VON HELM R LGB4 

15 ESTERO BY LWB 

8 ES LER CeCe war 

14 MORRO EY LBB1 

13. PT BUCHON LWB2 

10 “WESTDAHL fea 

S SOUZA R LGB 

12 LANSING BR B 

V1, VP SSAn EUs heen 

4 SANTA CRUZ MB 

2 AN ISLAND LWB8 

2 MONTARA LWBI1O0A 

1 YB ISLAND 
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G 


~1 Cy 


APPENDIX C 


LISTING OF BUOYTENDER ROUTE SELECTION PROGRAM 


PPewe eeROGRIM *** PUOYTENDER ROUTE SRLECTICG™ *** 


PeeircIT CHAPACTFR*15 (a) 
Mees, TDIS,TCRDIS,DIS (15) 
Pumeerr 1,J,L,C,2,NUM,NCD, CAM(15) 
CHARACTER*15 FNAME,OMNAME 7 
COMMON /NAME/ ANAMP(77) 
SOVMON /ValLkK/ N, 1TOUR(71) 
eeeiON /TOST/ TOC7S, 7%) 

Som ON /INVI/ iT (1 24) 

GOrMCN Joust, D7, 7) 
COMMON /SAMT/ ISTCR(1S,71) 
ao ON ie, TTPTOR(14¢) 


Sob OBTATN (FNAME, CNAME, D) 


OPEN(3,FILE=O™AME,STATUS='NEW' ,FORM='FORMATTED' ) 
TORDIS=99999,99 

CALL NODE(N,NUM,SAM) 

CALL NBR(NUM,SAM,DIS) 


Wee? L=1,5 


DOs a =) N+1. 
CeO )=aT STOR (ings) 
SOha Leue 
CALL MARFD() 
Di: ee a 
CALL TRPREP(C) 
R=( 
= 
Ph) =P POUR (1) 
in ( 2) FP ROur (2) 
Stes hee YY (i, 6, GoD) 
T=I[+l] 
Peek. ( COO) ) eGOTrol 1) 
CrhbeaADITOR (TT, RpuOb) 
GOT Ga. 
PEG 1 Ger ole 
CALL REVTOR(R) 
COT C 25 
ies) 6 
Dem ava) =i pN 
CPico=—Dle +p( bP POUR(J)e@ (TOUR(J+1)) 
COMTLN SE 
Pega, Ge. 1 OP DIS) GOT: / 
mere bs=i| DIS 
Power =i ‘+1 
Pome. (4 y=) TOUR (3) 
GUN a hae = 
aE, che 
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ae 


OY 


oO} ~ 


——m | 
rs 
C 


mm wh 


ee 


oy 


TP 
a 
AMRIT fat 
WRITE ( 
ick eee eral 
( 
( 
( 
( 
( 


wd 


Weil 
Has 


i Re 
ee 
CEC SF 
5S ber 
FOPMAT 
FORM ec. 
Pus kr 
ig Sasa 
Eh 


= % 


ITOUR (1) ,AMAME(T TOUR (ee) — tee 
(ITOUR (1) , SMAMEUT POU Ree — eo) 


= 


bh NO Ferm C2 SR 
eee eee eee eT 


(0 WW a Rad re ie aes 


~——— wy 


(1X,'ThHe SOlLUE16ON FCO= a ae eee oe, 
(1X,'Totil tour Aistates = ieee ee 
pe The tour is as ftollows as 
Cy op) ee) 


SURROUTINE CBTAIN(ENAME, OMAME, ‘!) 
IMPLICTID GyARiGTES* lore e 
INTEGER N 
REAL TIME 
CHAPACTER* 1] FNAME,OMAME 
So'MON /NAME/ SNAME (73) 

CIUMON /O1ST/ Begone 


a 

,7) ONA“E oa 
oe PORVM="E@e et er 
OU EE = eee 


eG) Se AOL ARN REDS. ARNARODAREN RRO 


ul aes 


’ 

TZ ((O (le eya iy) , l= ia 
7,CSR=L4) [Gee eta) 
) 


,'Untor waat €ile namo 1S eee oe 


ov 


Dc 
tel 
fez 
Ls) 


ho 


) 


fO 


Aa 


tO 
- 2 


FORMAT (13 
ceo Or ie ee, Wie eles aoc Oxide ely), 25.2,' sscon” 
PORMAT(LY, ‘to solve this 2robliem. Please wait.'///) 
See. (Pa. 2) 
meee i (ix, Frrvor ucon Bend” ) 
END 
SUS 2eUTINs NODE(N 8 U07,SA™) 
PEG ed, TTEMP,NUN., Sey je aie 
SES VON RAND? 1X;,1Y,12 
i= 1S 2G 
TY=2544% Ep fp, 
i201 79 
oun +h (ge * ALOG(ELCAT (t =4) )) 
CC ae eee = 1M 
ae ee er (ORNPOM( 1, Vee UGeAT (%.) 9 
Pole Poe. 9) GeO 3 
F=a 
DC 4 K= ie ne ie if ho f © vj 
LF(ITEMD. mee SCS LZ : 
ONT eR 
LEC pe peree TO? ee 
See a a ee 
SOR meena 3 
RETOUR. 
ae f 
ee@eett eae eo (NU, SAY, DIS) 
IMTECER 3, sy UMCLONS Sy vUsu zak Sip Gyr oii Malas renee 
Pe eee eo, OTS (15) 
Oe ae 8 K/ Ny, ITCUR(TT) 
Geese PI ST/ O(7C, 7%) 
Ge MOM eoAYT/ ISTOR(15,71) 
Se wien ey Pps e/ TO{(7C, 7%) 


eae Ga 
DP. ODR=S Al (Lb) 
Cimce l=), %. 


rr 
SA Ne 
L ee 


BO 2. =), %| 
Bleed} = 


Cx eee 0 
mac DE} 


2s rere 


ee ee 


= aN Ole Ts: 
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Ss ) 


Paes 


Ze 
en 


r+ 
a) 


NO to 


TEMPP=7D (12 ROEM et ein Bes 
De 72378. 
TE (TD (LETRA Tin) Cie te a 
17 chP=s 
TE4PD=Tb (lo Orie ae 
CONTINUE 
ISTOR(L, 2 Pleat Tetee 
DLS (Ll) =p iS te ae Gem 
DG 239) = 1 
TO(K, ITEME) =9909., 19 
SOT in, oe 
CONTINUE | 
LSTOR(L, ttl wet api. 1) 
Les (L)+D (ISTOP(L;N) 00S Se ee) 
CON TINGE 
Or x= 


CO. 2) ee i 
IF(SIS (1). FO. DIS (141) )01s(t ee 
IF (DIS(1). GP oe IS t1 +1) Cees 


NOEX="% 

TOMPD=D 1S (141) 

DIS(1 sly pee) 

Clot i= vee. 

BO. 22K J= been 
LTE“ P=1 STC rie io 
[IS TCR Ge bya So ee es) 
ISTGR(L,J)=ITE™P 


Ce ewe 
ON bade ee 
Pe (hOrx <P. 1)GOre 7 a 


> 


Rane Ro 
EWC 


SUBROUTINE TRPREP(C) 
Leas ot, 

COMMON /'s Hee © WITOUR (71) 
COMMON /TMTP/ [oon aera 
J=U 
WOtee) <= 17 

Loeler {har OUR) 
PTET nth) = 14 eee 
LF(ITPTOR(K).FO.C) J=K 
aren wmynnits 
DC 41 Ye] am 
PICCa (pal t eprce (ae) 
SO OG 


Nal 


oo 


oye 


() 
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TOUR (N+1)=1 TOUR (1) 
RFTLUPN 
END 


Seer euUlTITRE SoReP () 
ie? CER d,s 


SOMO) ole, < TTeus (7) 
Seo gost, TD(7*, 7%) 
DO 5% g=1,XN 


oi 
nen s 
Mm 
i 


UR Las 
Diy = — 1b a 
INU 


Poem Ueed), [TCUR( ITE) )=H=t 
Migempere (J+l), 1TOUR(I) )=1.% 
eee ec 
Pega wee) 7 fT CUR(]L))=1.9 
ener rai) pLTCUR(T))=Hlac 
Bae UN 


Ens DD 


SUBPOUTINE SE 
Pee eee to COL, 
Powe cs Pp fo TN ( 
aot, Cc, TEAR, SY 
COMMON /AALK/ ¥& 
een Test / T 
COE SS VN aay et 
Geer e/eist/ £ 


ieee smo, ft: 


Chae aaa 
CO 47 X=1,5 
Fa 
TEMP=9099,99 
DC 61 J=3,"! 
eo Se (2) remit) Ck 
Le (ieee (TT CURT Ieen Ch 
i=) (POs (2 ei Pour ia)) 
FoJ 
SOC Mars 
PeiGseee sO VGOTORA 
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a CO 


AY (1 TCUR CEs < 
Ke ee 

( 

( 


oe OO 


Ky =1T1e0rs) 
(K)=1 TOUR (F-1) 


KR) =PaL 


-U 


Ag a (Ae TAU) Cab ian 
[ee ee = 


ee) 


e, 
L? 


a 
Fi) 
— 
— 


NU aoe 
eo 


S 
ceil a. 
7 
a 
rx) 
as 
U 
| 
>) 


CN 

Lara See enee 
;CD=N0D (COL) 

IT (25+) H=S 7 ee) 
Poe lee2 aS PE) 
RETURN 

END 


SUBROUTINE ADGTOR (1 yay Ic 
[UPEGER. I ogee ice 
COMO JUNLX/ = rouse 
COMMe™./TDST/ 2Daan 
Cee /TUVT/ 1 ele 
COO /TMTE/ TPPren etc 
eXCo 


, 


ee 


ae wea 
eel, CRM [Obs 2=J)=1 teres 
COMTINCE 
Dm 5-7 168 
PPOUS(SY=1 Tea oom 
CONT INGE 
TOCLT (251 =) ) 5.0 oe ee 
TOUT (2 * 1=15 21324 ie) = 
TO Ceo * 1 = 1) Ce ee ee 
TOVLT(eF 1) 10 (tee 
Teer 2 #1.) To ene 
R=] 
PECUR. 
END 


Se, 


E)=D(ITOUR (1), 1TCUR (2 )3)- Oi eS 


OF 
Oy 
( 
A 
7M 
ves 

A(K)=SLA(K)4+D (1 TOUR (PF) GU Ge a ee 
( 
“ 
rt 
‘ai 
E 


Poe he Bes. 


) 


(i 


olak 


1, 


See OCT INE PEVTOE (FP) 
Pir GER FP 
eee ON PA ALK/ 17 f POUR (71) 
MeO ST TR TT ETC eR (14%) 
eS Ga), 41 
Pere ey = = OUP CI) 


ee Oe ee, 
Se. 12 
R=1 
he lor Ny 
ee, 


FUNCTION PA'DOM(A) 
Mee Ge Tee 

Bo Oy Veen’ IY, 1Y, tz 
eee eet ~, 1 77)-2* (TX/13 
fee ore OD (1, 176)-2* ( 

me os obit 717 8)-27* (17/178) 
IF(1X.LT.9) IX=1X+30269 


memory, LT eC) 1 Y=1 Y+2 0277 

IF(IZ. LT. CY IZ=TZ+30223 : 
TEM D= ee 2 *°. em Nm (i / 3027 7 
Belles en (TEMP, 1.9) 


ee 
orn 
a3 & HY UF 
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Tae ue (I 


fond 
“i 


Wi coe 


C) 


is2 


Is 6.3 NRO 


Rk te WD OO JON UI 


Ca 


AER Doe 


FeIPGIC PP Creve Ne eRe See 
PNTeECHR eras 

COMMON SE Ls vo er wee aa 
COMMONS INA IE (sie aa 
WET Ue Soper) 

He Re Lotd( aoe) 

ce Tea (se oes) 

i eee 

OS Alec een) 

RICE (ae) 

aN Pele eee espe) 

vc Po RBar ) 

S PolST Gee) 

SLT ET. Ce ee) 

oy Sens BE ala te ails eS 

IF( (TANS .2LT. 1) CR. (CIS NS. GT. oi ee ee 
Te (LANG Oe lk) CAL ele ae 
LE CERN Ss EQ, 2) CALL Pe ope ae 
ee es EOD, Cee SiS eC) 
EE (1ANS 22 2 4 OA eer Te) 


Pe (eee 
Pe Chases. 


POI S}CAtt 


APPENDIX D 


LISTING OF MATRIX UTILITIES PROGRAMS 


BNC E MAT Rs UO aig Ge ae 


iors &) 
VE sO) Ce Ole 


') 


STOP 

PORMAT(1X,'This is a sat of Gtilitves # ope ws-om ae 
ik) Ca aie 92rograr RCUTE. RUN a ney 30S fom eee 

7 OR ieee ee and mnanivoulate a Gistante=442¢n eee 

Oo a a, - create a t1ist3znce= natrix") 

PORFMAT(IM,'2 = extract a small Matrix fren a Yargene aa 
FORMAT(LX, "3 = Correchea Srey leu sao Nace ix 
EORMAT(IX,'4 - create a Rartcooy cl alae i7e ere 
FORM. T (1%, 'S - list names with their vow number in 7 ote 
PORWR. (LRy) 6 = See. in iS mci nae) 

SURMAT(//1X, “Nat 1S) you™r S]= Necro. a. 

FORMAT (I1) 

an B, 


“IMPLICT 
Rea i © 


Ree ae 


Cee Ie 
COMMON 


VATE Tee) 


IT CHARNC To si) 
Poo 

Pe 
JOS Te ae ae) 


JNA CoN Se (lee) 


* 
po) 
oe) 

ee oe) 
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) 


re. 


Sees | 


Ud 
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= OV Ol Go oa 


em UW WI LY WW WO W LW 


pe ENS 


a> 
IND 


©) 


Je 


DO 34 1 

WRITE (*, 
PFAD(*,2 

Gore tNuEe 


bi 


WRITF(*,37) “IN 
WRITE(*,32) 
Beeeo [=],%-] 
D(I,1I)=9999,99 
DO 40 J=I+1,N 
WPITE(*,41) I,J 
READ(*,42) DIST 
Dime =O1ST wa 
DOmIea L=1,100R 
COM RUF 


alr. ST 
PMi@eee te Gh. 0) GCTO42 
J=N 
I=N-1 
COT LNGE 
SOUT INUE 
D(N, Ni) =99°99,90 
CALL STOPE(N) 
RETURN 
PO Vwi oO. large 1S the problem? '\}) 
meRMAT (BN, f 2) 
Beet ieee Inout Names for noreport and aids') 
Boe lee Groen’, 13, ' Name? ‘\) 
FOPMAT (A) 
Sermo were ts," BY ',12,° Matrix’) 
POP ou Aista eee mao ne OSTme ho OO lnt J") 


Rebeewe(e mmo meStanc> from SO *,713,' to ne *,12,' = ‘'\) 
memes) (Bt,,67. 2) 
Se, 


SUBROUTINE EXTRCT() 
IMPLICIT CHARACTER*15 (A) 
INTEGER N,“,1I,7, [LROW(1LOS) 
REAL §D(107,1°%) 
CHARACTFR*15 ASMNM(1°7%) 
CGe er JDIST/ DN(109,1°%) 
CO™MOMN /NAME/ BNAME(L C7) 


Ve ee ye) 
reo belee(%5,.1 L ) 
Vie ee 1% ee 
amie. ( * (les: ) 
Cri OBT ees 
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eee ee) 
RFAD(*,15)™ 
DO 20 I=],* 
ee ei en 
eal FORMAT(1X%,12,'. ROw number Proumeen eee ee 
READ (7; 2S) ie ale, 
Per? = 
SOR Gy rele— 0 eC OR (ie Jt) 
oe COULINGE 
ASMNM(1)=ANAMF(IROW(I)) 
2h See wuss 
EO e222 J=t.k 
Ol l= 
D(I pes a ee oa ale 
Za SST Nae 


Dd a ee 
EA ia ea) 


oe eG aoe 
OPN in oie ees 
PETURN 
4 ECEMAT(//1%,'To create a snaller Matrix from your saan 
iP): PRM cr lee ee matrix, anter EF» 2Ow mtr eer sera 
see FORMAT (1K,'nain Yatrix aits wWhied 92 72367 6S pee 
I SCS MAC le eo ars gan epee . 
les POR MAT (jo ‘ Aat 18 the SizS Of the snaller sia... =a 
ES ne er 
END 
. 
Dee tiek <i ae ay 
[NM PETeIT CANR A ies on) 
Pet eGre 1, Jel eee eee 
Peto ia 


coMMON Panes 
SNAG See ey, 
ibe 2 et ee es) 
WRITE (*,18) 

16 Se Ss SUSEN SA ane) 79) 
VIPITE (*,12) 
eee) ee 


Oia eee Ce eo. 
1 es ae Sire a eee a ue es 
LON Sac se le Oe 


“RITE (*,21) 
SPS ON ey mete 
Risse a 


pee 
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Reve (.* 
WRITE ( 
PG ae 
SDs Ga 
ae) = 
BiG wees ace 
Be CONTINUE 
Swe, 1R)=Di1St 
Ge7O15 
fain (* 26) 
Pe ie (* ,27) 
So bi* 1218 
WRITE (*,23) 
Beam * ,)2) TC 
DC 33 [=IR,'N- 
Be 2ae7) =1C 7% 
apes, 2). OL dg { 1d) 
Veoh {* , 2 4 ) 
Pee (x 12) DIST 
Peat ser. i). COro 22 
TI =N 
J=N 
GC a: 
a Pee e- OST 
pomieage =), 1006 
1ia Sow eee 
Dio, 1 )=D1ST 
3] S265 a OES 
ee= T+ 2 
ea oC oO 


L 


a 


LPG Ot leer) 


Pe * 7 
™ 


ris 
— 


Moto po hoeryg 


~ 


DIST 
GS 
19%¢ 


ee) 
U1 


(J 
in 


poeta T (3h ,1 2) 

FORMAT(E7.2) 

Serer (//i<,°toi1s utility is for SS ane Saeco yoletila st +) 

Pee (EX, 2 5>revyiously Stores cata file 

Peri iy,’ | = Single sntries 2 

Pein (1X, ' Atk twcOrines bons Selec 

oe aL 1} 

Bot el, 

Bowe AT (le, 

Seo |< 

Coe Clef 
f 


a! 
SCOMAT (IN, "7. 


a ol el 
Voie ous fe Le eS a 
Np) 


t= 


I o> GiMeGe CG omieec bo’ \) 

>lumn HUNGER MemGOorrect * \) 

eee is aise con] hoy boo tay =, le 2) 
(eesca ace. = \) 

Bie Vet Gee sta bye 2 1Sc4t ive Sirsttacl” ) 
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ROUND ND ee 
Sy to tor 

2.0)8 PRS) 
oe 


NM ho 


C) 


es) 


coat © Ue iO 8 airy ee tor | 


faae 
CSS ta 


re? 


SUSSOURi rs 


FPURMAT (US, St age pemeroc 


ORG (le Se ape 


MATPRT () 
es CIN 
[Ne@gear 
COMMON 


Do Ge pee 


COM MOR air se 
YR cap Gd Ol ce 
iS, le, ell 
WRITE (*. 50) 
ELT Paes) 


CALL COTAIN (MN) 


OPEN (3,21 LES] At ee Tose 


SO] Uy A Tes te 


CHAR Ae ree ele 
age baie re a 
{EIST “Cee 
AL ee Lane) 


Ne) 


in he eee 


(A) 


u~ 


oN US a 


,»,FORM='?CRMAT Team 


TS=1 
TF=N 
i eae an = oo 
JS=] 
JE=9 
S/OWbce SOMA i 

IP (be. eee 

DOE g=gse = 

eA sei ie 

COM ae 
COE ee UP 
ies pe ee 
JS=J E+) 
JESHIF+9 
es ee ae J 
Leer 5 "ij Gane 
igs ont 
lle ee ata 
Le meer.) tae 
Porro) ace ayes 
Gtesn 4 
Vi ee oe et) 

RETURN 

Se eee utility Dra0areas a £11 >= 2onte ieee 
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